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Abstract 

Let s = (si, S2, ■ ■ ■ , Sjn) and t = {ti,t2, . . . , t^) be vectors of non-negative integers 
with Yl^i — Sj=i ^j- B{s, t) be the number of m x n matrices over {0, 1} 
with j-th row sum equal to Sj for 1 < j < m and k-th column sum equal to for 
1 < k < n. Equivalently, B{s,t) is the number of bipartite graphs with m vertices 
in one part with degrees given by s, and n vertices in the other part with degrees 
given by t. Most research on the asymptotics of B{s,t) has focused on the sparse 
case, where the best result is that of Greenhill, McKay and Wang (2006). In the 
case of dense matrices, the only precise result is for the case of equal row sums and 
equal column sums (Canfield and McKay, 2005). This paper extends the analytic 
methods used by the latter paper to the case where the row and column sums can 
vary within certain limits. Interestingly, the result can be expressed by the same 
formula which holds in the sparse case. 
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1 Introduction 



Let s = (si, $2, ■ ■ ■ , Sjn) and t = (t^, • • • , ^n) be vectors of positive integers with 
Y^iLi^i = 'YTj=i^j- Let B{s,t) be the number of m x n matrices over {0,1} with j- 
th row sum equal to Sj for I < j < m and k-th column sum equal to for 1 < k < n. 
Equivalently, B{s,t) is the number of labelled bipartite graphs with m vertices in one 
part of the bipartition with degrees given by s, and n vertices in the other part of the 
bipartition with degrees given by t. Let s be the average value of Si, S2, . . . , and let t 
be the average value of ti,t2, ■ ■ ■ ,tn. Define the density X = s/n = t/m, which is the 
fraction of entries in the matrix which equal 1. 

The asymptotic value of B{s, t) has been much studied, especially since the celebrated 
Gale-Ryser Theorem [10] that characterises such that B{s,t) > 0. Various authors 
have considered the semiregular case, where Sj = s for 1 < j < m and = t for 
1 < k < n. Write B{m, s; n, t) for B{s, t) in this case. For the sparse (low-A) semiregular 
case, the best result is by McKay and Wang [7] who gave an asymptotic expression for 
B{m, s;n,t) which holds when st = o(^{mnY^^) . In the dense (A not close to or 1) 
semiregular case, Canfield and McKay pLj used analytic methods to obtain an asymptotic 
expression for B{m, s;n,t) in two ranges: in the first, the matrix is relatively square and 
the density is not too close to or 1, while in the second, the matrix is much wider than 
high (or vice-versa) but the density is arbitrary. For the sparse irregular case, the best 
result is that of Greenhill, McKay and Wang [2], who gave an asymptotic expression for 
B(s,t) which holds when max{Sj} maxjtj^.} = o((Amn)^^^). 

See pp, [2] and [7] for a more extensive historical survey. 

The contribution of this paper is to adapt the approach of [1] to the dense irregular 
case when the matrix is relatively square and the density is not too close to or 1. See 
McKay and Wormald [8] for the corresponding calculation for symmetric matrices. 

In keeping with these earlier papers, the asymptotic value of B{s,t) can be expressed 
by a very nice formula involving binomial coefficients. We now state our theorem. 

Theorem 1. Let s = s{m, n) = {si, S2, ■ ■ ■ , s^) and t = t{m, n) = {ti,t2, . . . , t„) be vec- 
tors of positive integers such that YlJLi = Y2k=i for all m, n. Define s = Sj=i ^j' 
t = ■'T'"^ X]fc=i ^fc; X = s/n = t/m and A = |A(1 — A). For some e > 0, suppose that 
\sj — s\ = 0{n^^'^~^'^) uniformly for I < j < m, and \t^. — t\ = 0(m^^^^^) uniformly for 
1 < k < n. Define R = Yl^=ii^j "~ ^)'^ ^'^^ ^ ~ Y12=ii^k ~ t)'^ ■ Let a,b > be constants 
such that a + b < |. Suppose that m, n 00 with n = o(A^m^^^), m = o{A'^n^~^'^) and 

(1-2A)V 5m 5n\ 

H ^ < a log ra. 

8A \ 6n 6m J 
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Then, provided e > is small enough, we have 



- i'ZT ft (;) n (:) (4 (1 - ^) (1 - ^) . o(,r',) . 

j — 1 k — 1 

Proof. The proof of this theorem is the topic of the paper; here we will summarize the 
main phases and draw their conclusions together. The basic idea is to identify B{s, t) as a 
coefficient in a multivariable generating function and to extract that coefficient using the 
saddle-point method. In Section 2, equation ([T]), we write B{s,t) = P{s,t)I{s,t), where 
P{s, t) is a rational expression and I{s, t) is an integral in m + n complex dimensions. 
Both depend on the location of the saddle point, which is the solution of some nonlinear 
equations. Those equations are solved in Section [3l and this leads to the value of P{s, t) 
in fl20l) . In Section |H the integral /(s, t) is estimated in a small region 7^' defined in (|33|). 
The result is given by Theorem [2] together with fl2^ . Finally, in Section it is shown 
that the integral /(s, t) restricted to the exterior of TZ' is negligible. The present theorem 
thus follows from (P), ([20]), Theorems EHS] and □ 

Note that the error term in the above slightly improves the error term for the semireg- 
ular case proved in [1]. 

Thereom [T] has an instructive interpretation. Write it as B{s,t) = NP1P2E, where 

-G::)> ---'ft(::). ---'ft(;:). 

E = exp (--(1- ( 1 - + Oin-')] . 

^\ 2\ 2AmnJ\ 2AmnJ ^ ') 

Clearly, N is the number of m x n binary matrices with \mn ones. Pi is the probability 
that a matrix randomly chosen from this class has row sums s, while P2 is the probability 
of the similar event of having column sums t. If these two events were independent, we 
would have E = 1, E can be taken as a measure of their non-independence. For the 
case when s and t are vectors of constants, that is, i? = C = 0, Ordentlich and Roth [9] 
proved that E <1. 

It is proved in [2] that the same formula for B{s,t) modulo the error term also holds 
in the sparse case. Specifically, it holds with a different vanishing error term whenever 
max{sj}max{tk} = o{{Xmnf^^), R + C = 0{{Xmn)^^^) and RC = 0{{XmnY^^). In [T], 
evidence is presented that the formula is universal in the semiregular case {R = C = 0) 
and it is tempting to conjecture that the same is true in the irregular case for a wide 
range of R, C values. 

We will use a shorthand notation for summation over doubly subscripted variables. If 
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OC J ^ IS 3j variable for 1 < j < m and 1 < k < n, then 



n m m n 

EST" \^ \^ 

k=l j=l j=l k=l 

n—l m—l m—1 n—1 

•'^jki •^*fc / J jki / J / j jki 

k=l j=l j=l k=l 



for 1 < j < m and 1 < k < n. 

Throughout the paper, the asymptotic notation 0{f{m, n)) refers to the passage of m 
and n to oo. We also use a modified notation 0{f {171,11)), which is to be taken as a 
shorthand for O [f {m, n)n^^^^'^) . In this case it is important that the 0(1) factor is 

uniform over e provided e is small enough; for example we cannot write f{m,n)rS'^ 
as 0{f{m,n)) even though e"^ = 0(1) {e being defined as a constant). Under the 
assumptions of Theorem [T|, we have m = 0{n) and n = 0{m). We also have that 
8 < A'^ < 0{\ogn), so = 0(1). More generally, A'''m^^+^--''n''+^'^ = 0(n^2+^*) if 
Ci, C2, C3, C4, C5 are constants. 



2 Expressing the desired quantity as an integral 

In this section we express B{s,t) as a contour integral in (m + n)-dimensional complex 
space, then begin to estimate its value using the saddle-point method. 

Firstly, notice that B{s,t) is the coefficient of xl^ ■ ■ ■ Vi^ ■ ■ ■ ytr iii the function 

m n 

HHil + Xjy,). 
j=i k=i 

By Cauchy's coefficient theorem this equals 



where each integral is along a simple closed contour enclosing the origin anticlockwise. 
It will suffice to take each contour to be a circle; specifically, we will write 

Xj = QjC ^ and = r^e ^'^ 

for 1 < j < m and 1 < k < n. Also define 



for 1 < j < m and 1 < A; < n. Then 1 + x^y^ = (1 + qjrk){^ + \jk{(^'^^'^^''^ - 1)) 



so 



B{s,t) 



nr=inLi(i+g.^.) 
(27r)™+"nr=i^?nLi-t^ 

r r nr=inLi(i + A,.(e-(^-+^''^-i)) ^^^^ 



where 6 = {6i, . . . , 6^) and = . . . , Write i?(s, t) = P{s, t)I{s, t) where P(s, t) 
denotes the factor in front of the integral in ([T]) and I{s,t) denotes the integral. We will 
choose the radii qj, so that there is no linear term in the logarithm of the integrand of 
I{s,t) when expanded for small 6, 0. This gives the equation 



m n 



j=l k=l 3=1 k=l 

For this to hold for all 0, we require 

A,. = (1 < J < m), 
A.fc = tk {l<k< n). 

In Section [3] we show that ([2]) has a solution, and determine to sufficient accuracy the 
various functions of the radii, such as P{s,t), that we require. In Section H] we evaluate 
the integral I{s, t) within a certain region TZ defined in fl22|) . Section [5] contains the proof 
that the integral is concentrated within the region TZ. 



3 Locating the saddle-point 

In this section we solve ([2]) and derive some of the consequences of the solution. As with 
the whole paper, we work under the assumptions of Theorem [H 

Change variables to {aj^JLi, {bj^}^^i as follows: 

l + _ 1 + /ox 

qj — 2 ' "^k — 1"- 27~' ^ ) 

1 — r ttj 1 — r Ofc 

where 

Equation (j2]) is slightly underdetermined, which we will exploit to impose an addi- 
tional condition. If {qj}, {r^} satisfy ([2]) and c > is a constant, then {cqj}, {r^/c} also 
satisfy ([2]). From this we can see that, if there is a solution to ([2]) at all, there is one 
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for which XljLi < Ylk=i > 0, and also a solution for which ^Jli dj > and 
Ylt=i ^ follows from the Intermediate Value Theorem that there is a solution for 
which 

m n 

^"^(^j = ^'^bk, (4) 

j=i k=i 

so we will seek a common solution to ([2]) and (jlj). 
From ([3]) we find that 

AjfeZ X = l + aj + bk + Zjk, (5) 

where 

a^buil — r — r a,; — r hu) 

Zjk = 2 (6) 

1 + r ajbf. 

and that equations ([2]) can be rewritten as 

Si — s 1 . Z,- 



a. = J2h-— il<J<m) 



Xn n ^ n 

k=l 



1 z 

h = — {l<k<n) 

Am m ^-^ m 



t^ — t 



Summing ([7]) over j, A;, we find in both cases that 



va n 

'r>'^a'j + ^^h = -Z„ . 
i=i k=\ 



Equations (jlj) and (jHj) together imply that 



3=1 k=l 



Substituting back into (j7j), we obtain 

'^j = • • • 5 '^m5 • • • 5 ^n)) 

^fc = • • • , «m, • • ■ , K)^ 

for 1 < j < m, 1 < k < n, where 

S j 5 Z j ^ Z ^ 



. . . , a^, 6i, . . . , 6„) = + 

Efc(ai, . . . , a^, &i, . . . , &n) 



An n 2mn 
tk — t Z,j. Z,, 



(7) 



(9) 



Am m 2mn 
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Equation ([9]) suggests an iteration. Start with = llf^ = for all j, k, and, for each 
> 0, define 



where Z^, = Zj, = Zj,{a\ , . . . , , fe^ , . • • ^ ) ^ind similarly for Z.^ = Z,^ and Z„ = 
. We will show that this iteration converges to a solution of using a standard 
contraction- mapping argument. 

Recall that A^'^ = O(logn) under the assumptions of Theorem [1] (which we are adopt- 
ing throughout). This implies that = O(logn). Therefore, within the region A defined 
by I a,- 1, 1 6^ I < n^^^^ for all j, k, we have that 



„ = o[m ' ) and = o[n ' , 



which imply that, in the same region, we have 



Therefore, by the mean value theorem, we have for £ > 1 that 



max a'- — a,- + max oL oL' =o(m ' ) max a'- — a'- 

-|- o(?2~"'"^^) max \b^l^ — h^l 



provided {af"')} U {6^'"'^} U {af } U {6^'^} C ^. 
Applying the iteration once, we have 

af = (s, - s)/{Xn) and b^^^ = [t^ - t)/{Xm). 

Since {af^}, {bf^} and {af^}, {b^k^} lie inside ^A, we find by induction that {af^}, {bk^} 
lie in j^A for all Moreover, the iteration is Cauchy-convergent in the maximum norm, 
and the error in stopping at {af^}, {b^i^^} is at most maxj \af^ — af \ +maxfc Ife^^"* — fo^^""^"* |. 

When we carry out this iteration, we find that all the encountered af^ and b^^^ values 
are 0{n~^^'^). It helps to know that the following approximation holds in that case: 

Zjj^ = (1 — r'^)ajbk — r'^ttjbl — r'^a^b^ — r^(l — r^)a^bl + 0{n~^^'^). 
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Using the fact that Yl]Li <^f^ — Ylk=i ~ ^® ^^^^ 

Z« = -A<"X:(6<'y + 6(n-'), 



k=l 
m 



Therefore, 



J;^' . +d{n-^) (l<J<m), 

An A (1 - X)m n 



Am A (1 — A)m n 

Similarly, 

n n 

k=l k=l 

m m 

j=i i=i 

m n 

3=1 k=i 

which gives 

(3) ^ , js^ - s)C il-2\){s^-sfC 

Xn ^A^(l-A)mV A=^(l - A) W 
{1-2X)RC ^^^^^5/2^ (l<j<m). 



2A^(1 - A)^mV 



, (3) ^ tfc^ , (^A: - Q/? , (l-2A)(f,-f)^/? 

' Am ^A'(l-A)mV A^(l - A) W 

(l-2A)i?C ^5^^-5/2) (l<A;<n). 



(11) 



2A^(1 - X)Wn^ 

Further iterations make no change to this accuracy, so we have that a,- = af^ +0{n~^^'^) 
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and 6fc = b^j^^ + 0{n ^^'^). We also have that 

(1 - 2A)(s, - - t) {s^ - s){t, - tf (s, - sf{% - t) 



Z 



\^{l-\)mn \^{l-X)m^n \^{l-\)mn^ 



{l-2X){s^-s)\h-tY ^i-2X){s^-s){h-t)R 

(l-2A)(.,-.)fa-t)C ~ 5/,^ 
^ A^(l-AfmV + ^ 



A sufficient approximation of Xj^ is given by substituting (fTTl) and f|T2l) into ([5]) 
In evaluating the integral I{s,t), the following approximations will be required: 



, (l-6A + 6A^)(g, -g)fa-t) , _3/2 
A(l-A)mn '^^'^ ^' 

(l-6A + 6A^)(sj- -s) 



2 



A,fc(l - A,fc)(l - 2A,fc) = A(l - A)(l - 2A) + 



n 

i2\ 



(l_6A + 6A^)fa-t) ^ ~ 1, 



Then 



j=i fc=i 



m n f ra n \— 1 

j=lk=l ^j=l k=l 



m n 

t 



j=l k=l j=l k=l 



^k' 



using ([2]). Therefore the factor P{s,t) in front of the integral in ([T]) is given by 

P{s,t) = (27r)^("+")A-\ 



(13) 



(14) 



A,,(l - A,,)(l - 6A,fc + 6A;,) = A(l - A)(l - 6A + 6A^) + 0{n-'^'). (15) 

3.1 Estimating the factor P{s,t) 

Let 

m n 
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We proceed to estimate A. Writing A^;, = A(l + Xjj^), we have 



(16) 

I ^ ^2 A(1-2A) 3 I A(1-3A + 3A^) 4 , ^( 4 ^ 
^2(1-A) ^'^ 6(1 -A)^ ^' 12(1 -A)^ ^'"^ V(l-A)' 

We know from ([2]) that A,, = mnX, which imphes that x,, = 0, hence the first term on the 
right side of ffTB]) does not contribute to A. Now using we can write Xjj^ = aj + b,^ + Zj^ 
and apply the estimates in (ITT!) and (fT2l) to obtain 

A = (A\l - A)^-^r"expf + + ~ _ (1_Z^^3 

^ ' "^KAAn AAm %J^rnn lAJ^n 2AA^7n 



(17) 

(l-3A + 3A^)i?4 , (1-3A + 3A^)C4 , ,1/2^ 
96A^n^ 96A^m^ 



where i?^ = X^jlil-^j ~ ^'^d = X]fc=i(^ ~ '^kY for any £. Note that R2 = R and 

To match the formula from the sparse case solved in [2J, we will write f[T7|) in terms 
of binomial coefficients. First, by Stirling's expansion of the logarithm of the gamma 
function, we have that 

TIT \ / x+d/-i \l—x—d\—N 

N \ [x (1 — x] 



X + d)Nj 



1-2X d^N {l-2x)d {l-AX)d^ 



{l-2x)d^N (1-6X)/X ^(d^N d 1 



24X2 gg^3 V X^ X^X X^X^ 

as X — > 00, provided x = x(X), X = X(X) = |a;(l — x) and d = d{N) are such that 
0<x<l, 0<x + (i<l and provided that the error term in the above is o(l). From 
this we infer that 

TT / ^ 1 TT / 1 (A (1 — A) ) 



Amn; J-^-L \SjJ \tj (4^^)(™+n-l)/2^(n-l)/2^(m-l)/2 



1 ^ J/ k=l 

xexp(-^-^-l^r!^ + ^U^^r- + -^ (19) 
' 4An AAm 2AA \n m) IQA^ \n 

l-2A/i?3 C3\ l^^fR^ ^C,^ ~ 

24A V n m / 96^4 m / 
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Putting (fTTI) and (fT9|) together, we find that 
P(s,t) = A-^(27r)-(™+") 

l(m+n-l)/2 (n-l)/2 (m-l)/2 / \-l m / s n 



2^(m+n+l)/2 y ii y^.J ii J (20) 

,l-2Afm n\ RC I - AA / R C\ ^,-in,\ 



4 Evaluating the integral 

Our next task is to evaluate the integral J(s, t) given by 

It is convenient to think of 6j, (p^ as points on the unit circle. We wish to define 
"averages" of the angles 6'^, 0^. To do this cleanly we make the following definitions, as 
in [1]. Let C be the ring of real numbers modulo 27r, which we can interpret as points 
on a circle in the usual way. Let z be the canonical mapping from C to the real interval 
(— TT, 7r]. An open half-circle is = (t — 7r/2, t + 7r/2) C C for some t. Now define 

= \^ = {x = {xi, . . . , Xj^) G I Xi, . . . , a;jv G for some t G M }. 

t 

If X = {xi, . . . , Xjsf) G Cq then define 

■ iV 

-y 



x = z-\ ^2^Z{Xj, 



More generally, if a; G then define x = t + [xi — t, . . . ,Xj^ — t). The function x ^ x 
is well-defined and continuous for x G C*^. 

Let 71 denote the set of vector pairs {d,<p) G C"" x C"' such that 

\e + ^\< (mr^)-'/'+2^ 

l^^-l < (l<j<m), (22) 

|4l < (1 < A; < n), 

where 9j = 9j — 6 and 0^ = 0^ — 0. In this definition, values are considered in C. The 
constant e is the sufficiently-small value required by Theorem [TJ 
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Let I^"{s,t) denote the integral I{s,t) restricted to any region TZ" . In this section, 
we estimate I^'{s,t) in a certain region TZ' 3 TZ. In Section O we will show that the 
remaining parts of /(s, t) are negligible. We begin by analysing the integrand in TZ, but 
for future use when we expand the region to 7Z' (to be defined in fl33|) ). note that all 
the approximations we establish for the integrand in TZ also hold in the superset of TZ' 
defined by 

\e + ^\< 3(mn)"^/^+2", 

l^jl < 3n"^/^+^ (l<j<m-l), 

\L\ < 2n-'/'^'\ (23) 
|4|<3m"^/2+^ (l<A;<n-l), 
I0J < 2m-^/^+^^ 

Define 6 = {§1, . . . ,9j^_i) and = (0^^, . . . , 0„„i). Let be the transformation 
Ti(^, 0, u, 6) = {0, 0) defined by 

u = e + cj), 6 = 0- (j), 

together with 9j = 9j — 6 {1 < j < m — 1) and 0^ = 0^ — (1 < A; < n — 1). We also 
define the 1-many transformation Ti by 

5 

After applying the transformation to Ifi{s, t), the new integrand is easily seen to be 
independent of 6, so we can multiply by the range of 6 and remove it as an independent 
variable. Therefore, we can continue with an (m + n — l)-dimensional integral over S such 
that 7^ = T;(5). More generally, if 5" C (-lyr, ^tt)"'^''-^ x (-27r,27r] and TZ" = T^{S"), 
we have 

/^"(s, t) = 2TTmn / G{e, 0, z/) dOd^du, (24) 

where G'(6>, 0, u) = F{T^{0, 0, i/, 0)) with F{0, 0) being the integrand of (pTj) . The factor 
27117171 combines the range of 6, which is An, and the Jacobian of T^, which is mn/2. 

Note that S is defined by the same inequalities (122!) as define 7^. The first inequality 
is now \v\ < imn)^^^'^^'^^ and the bounds on 



m— 1 n— 1 

0fc 



i=l k=l 



still apply even though these are no longer variables of integration. 
Our main result in this section is the following. 
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Theorem 2. Under the conditions of Theorem\^ there is a region S' ^ S such that 

1 l-2A/m n\ 1/1 1\ /i? C 



X exp — — — + — +TT — + - - + 

1-8A/R C 



IbA m ^ / 



In the region 5, the integrand of fl2^ can be expanded as 



- Yl + "^^kK^ + + ~ ' Yl 5Z(^3 + f^i'''^^^ + + 

j=l k=l j=l k=l 

m n m n \ 

j=l k=l j=l k=l ^ 

Here a^^, and 7^^ are defined by 

^^jki^ ~ ^jk) = A + a^-fc, 
iA,fc(l - A,,)(l - 2A,fc) =A, + (25) 

where 

A = iA(l - A), A3 = iA(l - A)(l - 2A), and A^ = ^X{1 - A)(l - 6A + 6A^). 
Approximations for ajj^, Pjj^, 'jjj^ were given in f[T3l) - f|T5l) . 

4.1 Another change of variables 

We now make a second change of variables {0, (p, z/) — ^2(C) '^)? where ^ — ((^i, . . . , 
and ^ = (^1, . . . ,^n-i); whose purpose is to almost diagonalize the quadratic part of G. 
The diagonalization will be completed in the next subsection. The transformation T2 is 
defined as follows. For 1 < j < m — 1 and 1 < k < n — 1 let 

= (j + ^k = ^k + dpi, 

where ^ 

and d 



I 1/2 I 1/2 
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and, for 1 < h < A, 

m—l n—1 

TT/i = Cjy Ph = Cfc ■ 

j=l k=l 

The Jacobian of the transformation is (mn)^^^'^. In [Ij, this transformation was seen 
to exactly diagonahze the quadratic part of the integrand in the semiregular case. In 
the present irregular case, the diagonalization is no longer exact but still provides useful 
progress. 

By summing the equations 9j = (j + ctti and (pk = + dpi^ we find that 



m— 1 

1/2 n 11^1/2 -1/2+e 



n-l 



(26) 



k=l 



where the right sides come from the bounds on 6'^ and 0„. This implies that 

C. = ej + d{n-') (l<J<m-l), 
^, = ^, + d{n-') {l<k<n-l). 



The transformed region of integration is T2 ^{S), but for convenience we will expand it a 
little to be the region defined by the inequalities 

|g|< 1^-1/2+. {l<j<m-l), 

Kfcl < (l<A;<n-l), 

I I ^ 1/2 -1/2+E /oc^ 
\tti\ < m ' n ' , [26) 

I I ^ 1/2 -1/2+e 
\Pi\ < n ' m ' , 

11^/ \-l/2+2e 

< [mn) ' 

We now consider the new integrand Ei = exp(Li) = G0T2. As in [JLj, the semiregular 
parts of the integrand (those not involving a^^, Pjf. or 7^^) transform to 

— Amnu'^ — AmT2 — Amp2 — 3iA^ni''K2 — 3iA^mi'p2 + 6A^tt2P2 

~ -1/2 (-^^^ 

— iA-^rni^ — iA^np^ — 3iA^cmTi7!'2 — 3iA-^dmpiP2 + A^nir^ + A^mp^ + 0{n ' ). 

To see the effect of the transformation on the irregular parts of the integrand, write 
Cm = ~ CTTi and in = — dpi- From ( l26|) we can see that Cm = 0{n~^^'^) and 
^„ = 0(n~^^^). Thus we have, for all 1 < j < m and 1 < A; < n, C + ^a: = 0{n~^^'^) and 
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CTCi + dpi + V = 0{n Recalling also that a^fc, /^j^, = 0{n ^^'^), we have 

m n 

^^ajkij^ + Gj + fi^kf 

i=l k=l -m n 

j=l k=l 

m n m n 

Y.Y.f^^k^'' + ^3 + ^kf = E + ^^)' + o{n-'/'), 

j=l k=l j=l k=l 

m n 

j=i k=i 

Moreover, the terms on the right sides of the above that involve Cm or C,n contribute only 
0{n~^^'^) in total, so we can drop them. Combining this with fl29l) . we have 

Li = —Amnv^ — Am:2 — Amp2 — 2tiA-^nv'K2 — 3iA-^mup2 + QA^tt2P2 

— iA^uTT^ — iA^np^ — 3iA^cmTiTT2 — 3iA^d'mpip2 + A^nit^^ + A^mp^ 

m—l n—1 

+ ^fe)' + 2(Cj + ik){i^ + CTii + dp,)) (30) 

3=1 k=i 

m— 1 n—1 

j=i k=i 



4.2 Completing the diagonalization 

The quadratic form in is the following function of the m + n — 1 variables ^, J^- 

Q = —Amnv'^ — AmT2 — Amp2 

m—l n—1 /Q-| N 

j=i k=i 

We will make a third change of variables, {C,^,^^) = ^al'^? ''") A*)) that diagonalizes this 
quadratic form, where cr = {a,, . . . , <7^_i) and r = (ti, . . . , r„_i). This is achieved using 
a slight extension of [51 Lemma 3.2]. 

Lemma 1. Let X andY be square matrices of the same order, such that exists and 
all the eigenvalues of X'^Y are less than 1 in absolute value. Then 

(/ + YX'^y^^^ {X + Y){i + x'^Yy^^^ = X, 

where the fractional powers are defined by the binomial expansion. □ 
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Note that and FX"^ have the same eigenvalues, so the eigenvalue condition on 

X^^Y applies equally to YX^^. If we also have that both X and Y are symmetric, then 

E (7) (y^-r = E (7) i^-yy' = E (7)^-'^ 

r>0 ^ ^ r>0 ^ ^ r>0 ^ ^ 

SO (/ + yX~^)~^^^ is the transpose of (/ + X^^Y)^^^"^ . Let V be the symmetric matrix 
associated with the quadratic form Q. Write V = + where has all off-diagonal 
entries equal to zero and matches V on the diagonal entries, and V-^^ has all diagonal entries 
zero and matches V on the off-diagonal entries. We will apply Lemma [1] with X = 
and Y = V^^. Note that is invertible and that both and V^j are symmetric. Let 
T3 be the transformation given by T^{cr, r, fi)^ = (C, = (-^ + KT^Kid)"^''^!''') ''"^ A*)"^- 

If the eigenvalue condition of Lemma [1] is satisfied then this transformation diagonalizes 
the quadratic form Q, keeping the diagonal entries unchanged. 

From the formula for Q we extract the following coefficients, which tell us the diagonal 
and off-diagonal entries of V: 



[C|] Q = 


—An — (1 + 2c)aj^, 


[il] Q = 


—Am — (1 + 2d)a^k, 


[A Q = 


—Amn, 








—2aji^ — 2daj^ — 2ca^fc, 






[C,^] Q = 


2o!j^, 


MQ = 





Using these equations we find that all off-diagonal entries of V^j are 0{n~ ' ), except 
for the column corresponding to u which has off-diagonal entries of size 0{n~^^'^). Simi- 
larly, the off-diagonal entries of KidKT^ all 0(n~^^^), except for the row corresponding 
to z/, which has off-diagonal entries of size see that these conditions imply 

that the eigenvalues of KT^Kid ^ire less than one, recall that the value of any matrix norm 
is greater than or equal to the greatest absolute value of an eigenvalue. The 00-norm 
(maximum row sum of absolute values) of KT^Kij is 0{n~^^'^), so the eigenvalues are 
all 0(n-^/2). 

We also need to know the Jacobian of the transformation T3. 

Lemma 2. Let M be a matrix of order 0{m + n) with all eigenvalues uniformly 0{n~^^'^). 
Then 

det(/ + M) = exp(trM- |trM^ + 0(n"^/^)). 
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Proof. The eigenvalue condition ensures that the Taylor series for log(/ + M) converges 
and that 

det(J + M) = exp(trlog(/ + M)). 

Expanding the logarithm and noting that ItriVf^l = 0{n^^^~'^^^'^) for r > 3 gives the 
result. □ 

Let M = V^^^Kid- As noted before, the eigenvalues of M are all 0{n^^^'^) so Lemma[2] 
applies. Noting that tr(M) = and calculating that tr(M^) = 0{n~^), we conclude that 
the Jacobian of T3 is 

det((/ + M)~^/2) = (det(/ + M))"'/' = 1 + d{n-'^^). 
To derive T3 explicitly, we can expand (/ + Kf^Kid)"^^^ while noting that a 



J* 



0{'n'^'^~^^) for all j, a^.^ = 0{m^^'^^^) for all k, a^^ = 0{mn^^ + nm^'^), R < mn^'^'^^ and 
C < nm . 



This gives 



i =1 



fc=i 



(^-^'5(„-')).H-0(„-), 



i=i ' k'=i 
+ 

m—l n—l 

j=i k=i 

for l<j<m — l,l</i:<n — 1. 

The transformation T^^ perturbs the region of integration in an irregular fashion that 
we must bound. From the explicit form of T3 above, we have 

m— 1 n—l 

= C, + E Oin-'/%^ + ^ Oin-'^')^, + 0[n-'l^)v + 0[n-^) = C, + 0(n-^), 

m—l n—l 

i=i fe'=i 
for 1 < j < m — 1, 1 < A; < n — 1, so cr, T are only slightly different from ^, ^. 
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For /i versus v we have 

where the second step requires our assumptions m = o(y4^n^^^) and n = oi^J^m"^'^'). This 
shows that the bound \v\ < (mn)~^^^^^^ is adequately covered by < 2(mn)^^^^^^^. 

For 1 < h < 4, define 

m—l 1 

Eh \^ /i 

j=l fc=l 

From (128|) . we see that Itt^I < m"^^^n~^^^~*'^ and |pi| < m~^^'^~^^n^^'^ are the remaining 
constraints that define the region of integration. We next apply these constraints to 
bound fii and z/^. From the explicit form of T3, we have 



. 2An 

m—l n—l , , , m—l 



i=i k=l j=l 

ca^:^, -1/2+e . doi** -1/2+e 1/2 I '^** 

n—l / ^ \ m—l 

fc=l /=i 

+ (32) 

= vTi + 0{A^^mn'^^^') 

I / 1/2 -1/2+5£/2n 

= 111 + o[m ' n ' ' ). 

To derive the above we have used 1 + c(m — 1) = m^^'^ and the bounds we have established 
on the various variables. For the last step, we need the assumption m = o{A^n^^'^), which 
implies that A~^mn~^~^'^^ = o{m^^'^n~^^'^^^^^'^\ 

Since our region of integration has \txi\ < m^^'^rr^^'^^^ , we see that this implies the 
bound l/iil < m^^'^rr^^'^'^^^ . By a parallel argument, we have 

z.i = p, + o(m-^/^+^^/V/2), 

which implies < n^^^m^^^^^^^. Putting together all the bounds we have derived, we 
see that 
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where 

Q = { \a^\ < 2n"^/^+',j = l,...,m - 1} n { iTfcl < 2m"^/^+', A; = 1} 
n{|/i| <2(mn)~^/2+'^}, 
M = { l^il < m'^'n-'/''-'' }n{\u,\< n'/'m-'^'"-''}. 

Now define 

S' = T,{T,{QnM)), 

n' = T^{s'). ^ ' 

We have proved that S' 3 S, so it is vahd to take S' to be the region required by 
Theorem [21 Also notice that TZ' is contained in the region defined by the inequahties fl23p . 
As we forecast at that time, our estimates of the integrand have been vahd inside this 
expanded region. It remains to apply the transformation T^^ to the integrand fl30|) so 
that we have it in terms of {cr,T,fi). The explicit form of T^^ is similar to the explicit 
form for T3, namely: 



j' = l k=l 



+ 0(n-^))/i + 0(n-'), 



\2An 

m— 1 7 , n— 1 



i=i fc'=i 

m—l n—l 
j=l fc=l 

for 1 < j < m— 1, 1 < k < n — l. In addition to the relationships between the old and new 
variables that we proved before, we can note that = /i2 + 0(r2~^''^), p2 = 1^2 + 0{n~^^'^), 
TTg = + 0{n~^), P3 = 1^3 + 0{n~^), = /i4 + 0(n~^''^), and P4 = 1^4 + 0{n~^^'^). 

The quadratic part of Li, which we called Q in ( l3Ti) . loses its off-diagonal parts ac- 
cording to our design of T3. Thus, what remains is 

m— 1 n—l 

—Amnp' — ^^(y4n + (1 + 2c)aj^)cr| — ^^(Am + (1 + 2(i)a;^fc)r^ 
j=i k=i 

m—l n—l 

= —Amnfi^ — AnfX2 — Amv2 — — Oi*kTk + 

i=l A;=l 
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Next consider the cubic terms of L^. These are 

— 3iA^nh'7C2 — 3iA^mi'p2 — iA^mr^ — iA^np^ 

m— 1 n— 1 

j=i k=i 

We calculate the following in Qn Ai: 

■ m—l n—l \ 



-'iiA.^nvTx^ = -?>iA^npp2 + ^^4^ ( + ^ ^*k^k j + 0(r 

^ j=l k=l ^ 

i.i'=i 

m— 1 n—l N 

j=i fc=i ^ 

Q • /I 2 "i-l 

3^ Trip2 \ — ^ — 1/2\ / ^\ 

-3iA3cmTi7i2 = -3iA^cnpiP2 H 2^ "j^o^i + 0[n ' ), (34) 

m— 1 n—l m—l n—l 

j=l A;=l i=l fc=l 



and the remaining cubic terms are each parallel to one of those. The proof of flMj) is 
similar to the proof of fl32l) . 

Finally we come to the quartic part of which is 



QA^'n2P2 + A^n-n^^ + A^mp^ = QAiP2^2 + A^np^ + A^mv^ + 0(n 



-1/2n 



In summary, the value of the integrand for (cr, r, /i) G Q fl is exp {L2 + 0{n ^^'^)) , 
where 

m—l n—l 

L2 = —Amnp — Anp2 ~ Ami'2 — 2, ~ / , '^*k'^k + QA^p2^2 

j=i k=i 

+ A^np4^ + A^mi'4^ — iA^np^ — iA^mi'^ — 3iA^cnpiP2 — 3iA^dmi>ii>2 

m—l n—l (Ofi\ 

- 3iA^npp2 - 3iA^mpv2 - i E ~ ^ E ^^-^^^ 

j=i k=i 

m—l n—l m—l n—l 

+ i E ^Ji'^J^/ + * E Kk'^-^l' + ^ E + Vjk(^]n), 

j,j'=l k,k'=l i=l fc=l 
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with 



no 



"^■((1 + cm + c'm^)aj^, + cma-'^) = 0{n ^''^^^), 



2Am 



Kk' = + dn + (fn)a^k + dna^^^) = 0{m ^/^+'), 

= ^ ^ {"majk + (1 + cm)a^:^. + dmaj^) — 3/3^-^ = 0(m^^^^^^^ + n^^^^^^^) 
Note that the 0( ) estimates in the last four hnes are uniform over k, k'. 



4.3 Estimating the main part of the integral 

Define E2 = exp(L2). We have shown that the value of the integrand in Q fl is 
El = £2(1 + 0{n^^^^)Y Denote the complement of the region A4 by Ai'^. We can 
approximate our integral as follows: 



Ei= E2 + 0{n-'/^) / \E, 

QnM JQnM JQnM 
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2 

QnM 



E, + 0(n-'/^) I \E, 



= [ E2 + 0{1) [ \E2\+0{n-'/') ! \E2\. (37) 

It suffices to estimate the value of each integral in fl37j) . 

We first compute the integral of E2 over Q. We proceed in three stages, starting with 
integration with respect to fi. For the latter, we can use the formula 



I x-l/2 + 2e 

ymn) 

■ (3' 



j expl^—Amnfi^ — iPfijdfi = j exp^- 



AAmn 



-ymn) 



provided /5 = o{A{mn)^^'^'^'^^\ In our case, /3 = 3^43 (n/i2 + ^1^2), which is small enough 
because of the assumptions m = o{A^n^~^'^) and n = o{A'^m^^^). Therefore, integration 
over n contributes 

/ 71 \i/2 / -QAl{nfi2 + mu2Y , ^/ -1- 



\Amn/ ^\ AAmn ~^ ^ / ^ 
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The second step is to integrate with respect to cr the integrand 

exp(- Anij,2 - ^ aj*cr| - ^^/^2 - i^nfi^ - 'iiA^cn^ii^i^ 

m—l m—1 m—1 n—l 

- ' + ' 3rj'<^j<y} + ^YY iVi'^l + ^'jk<y]rk) (39) 

3=1 jj'^i j=l k=l 

l2. 



This is accomphshed by an appeal to Theorem SI presented in the Appendix. In the 
terminology of that theorem, we have N = m — 1, 6{N) = 0{n~^), e = ^e, e" = ^e, 
e" = Se, e = 6£, and ^(A^) = £ + o(l) is defined by 2vr^^'^^'^ = N~^^'^^^. Furthermore, 

2 An ^ ( ^ ^ 

A = — — -, ttj = -aj^ + ( 6A4 - — — )i^2 + i 2^ VjkTk, 



m-V ' V ^ 2A 

^ _ An p _ QAln 

m — 1 ' AAm ' 

n — 1 

Jj = i Uj^.T^ . 
k=l 

We can take Z\ = |, and calculate that 



k=l 



Ei^ 1 \ ^ ni f 3A4 9^3 \ ~, „i. 



^ 1 ^ 1 1 

^ E ^j + E ^? = + E("i*)' 

j=i j=i j=i ^4Q) 



m / 3A4 9A^\ i pi, -i/2\ 

n \ A AA^J 2An " ^ ^ 



fc=i 

Z = Z, = e.p^-^+0{n )J=0(l)exp(^^^^J. 
Applying Theorem |U we see that 6*2 = 0(n^^^^), and so integration with respect to 
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contributes a r-free factor 



/ TT \(™-i)/2 /m/3.44 15Al\ 1 

+ 77^ E + ^(^"'^') + 0{n-'/'Z,) . 
4A n ^ J 

By the conditions of Theorem [H < n^^^ so 0(n^^^^) + 0(n~'^^^Zi) = 0(n"^''^) = o(l) 
as required by Theorem HI 

Finally, we need to integrate over r. Collecting the remaining terms from (!36l) . and 
the terms involving r from (!38l) and (l40l) . we have an integrand equal to 

n— 1 n— 1 . n— 1 n— 1 \ 

A;=l fc=l fe=l k,k'=l 

In the terminology of Theorem HJ N = n — 1^ ^{^) = 0{n~^^'^), e = |e, e" = |e, 
e'" = 3e, e = 4e, and i{N) = e + o(l) is defined by 2m^^^'^^^ = N~^^'^^^ . Furthermore, 



Am „ ?)A^m QA^m 

n-l An AA n 

P _ A^m p _ QAlm 



We can take ^ = f again and calculate that 



l^^^ + TJ^2l^Fkk' = -[-r^-TrTz]+0{n ] 



«=! fc,fc =1 

15^., 3^^^ 1 3A^,n 



«=1 A;,fc =1 k,k,k =1 

1 ^ 1 ^ 1 ^ n-l 

8^4^ 2A^ 

Z = Z, = exp^^ + 0(n ) J = 0(1) exp ^^^^ 
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We again find that 02 = 0{n ^^'^). Including the contributions from (138|) and (HTjl 
we obtain 



^ \Amn/ \AnJ \AmJ 



-l)/2 



3^4 /m ra\/3A4 15 
X expl 1 + — t + { — + —)( — I 1 

11 1 



(42) 



2Am ' 2An) ** AA^m 

k=l 

7 = 1 ^ 



Using f[T^ and the conditions of Theorem [T], we calculate that 

2 \n m/ 



m— 1 



Y,ic^^.f = l{l-2XfR + 0{n 



3/2) 



n-1 



3/2^^ 



k=l 



Substituting these values into fH2]) together with the actual values of A, A^, A^^ we con- 
clude that 



C ^ / TT \l/2/ TT \(m-l)/2/ yr \(n-l)/2 



,1 l-2A/m n\ 1/1 l\/i? C\ 
^^Pi 2 24A Vn m) 4AVm n/Vn m/ ^ ' 
\-%A(R C\ 



16 A^ 777,^ 

We next infer a estimate of Jq\E2\. The calculation that lead to P2|) remains valid if 
we set all the values ^3, Pjf,, g--i, h/^f^', Uj^ and Vji. to zero, which is the same as replacing 
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L2 by its real part. Since |-E'2| = exp (Re(L2)), this gives 



15Ai (m 



^2 = exp — ^ + % — + — + oil) / ^2 



= / E2 (44) 

JQ 

under the assumptions of Theorem [H The third term of (1371) can now be identified: 

d{n-^l^) [ \E,\ = d{n''/')n^ [ E, = 0(n-') / E„ (45) 
Jq Jq Jq 

where, as always, we suppose that e is sufficiently small. 

Finally, we consider the second term of ( |37ll . namely 



\E. 



2h 



which we will bound as a fraction of |i?2| using a statistical technique. The following 
is a well-known result of Hoeff ding ^ . 

Lemma 3. Let Xi,X2, . . . he independent random variables such that EX^ = and 
< M for all i. Then, for any t > 0, 

/ ^2 



Prob(^^Xi > < expf- 

i=l ^ 



Now consider IE2I = exp(Re(L2)). Write = Mi fl A<2, where Mi = { \fii\ < 
m^^'^rT^^'^^'^'^ } and M2 = { I'^il ^ n^^^m^"*^^^^^^}. For fixed values of [i and cr, Re(L2) sep- 
arates over T]^, r2, . . . , r„_i and therefore, apart from normalization, it is the joint density 
of independent random variables X^, X2, . . . , which satisfy EX^ = (by symmetry) 
and |X/j| < 2m~^^^'^^ (by the definition of Q). By Lemma [3l the fraction of the integral 
over T (for fixed /x, cr) that has > ri'^'^mT^^'^^'^'^ is at most exp(— m^^/2). By symmetry, 
the same bound holds for Vi < — n^^^m"^''^'^^^. Since these bounds are independent of fx 
and cr, we have 



/ IE2I < 2exp(-m^72) / {E^]. 
JQnMi JO. 



By the same argument. 



IE2I < 2exp(-n 72) / IE2I. 
QnMi Jq 
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Therefore we have in total that 



/ IE2I < 2(exp(-m^72) +exp(-n^72)) / IE2I 

< 0{n~') [ E,, (46) 

JQ 

as for (05]). Applying ([37D with (03]), (gS]) and (geD, we find that is given by fH3|) . 

Multiplying by the Jacobians of the transformations T2 and T3, we find that Theorem [2] 
is proved for S' given by (133|) . 



5 Bounding the remainder of the integral 

In the previous section, we estimated the value of the integral I-ji'is, t), which is the same 
as I{s,t) except that it is restricted to a certain region IZ' ^ IZ (see fl2T| - [23|) ) In this 
section, we extend this to an estimate of I{s,t) by showing that the remainder of the 
region of integration contributes negligibly. 

Precisely, we show the following. 

Theorem 3. Let F{6,<j)) be the integrand of I{s,t) as defined in (1^ . Then, under the 
conditions of TheoremUl 

I \F{e,(t))\dOd(t) = 0{n-^) [ F{e,(f>)d6d(f>. 

Jn" Jn' 

For 1 < j < m, 1 < < n, let Ajj^ = A + a^^ = |Ajfc(l — Xjj^) (recall (125!) ). and define 
Ajnin = miUj-fc Aji^ = A + 0{n~^^'^). We begin with two technical lemmas whose proofs are 
omitted. 

Lemma 4. 

m n 

\F{d,ct))\ = \[\[f^,{e^ + cp,), 
j=i k=i 

where 

fjk{z) = a/1 - 4Ajfc(l - cosz) . 

Moreover, for all real z, 

< fjk{z) < exp{-Ajkz'^ + j^AjkZ^). □ 
Lemma 5. For all c > 0, 

/87r/75 
exp(c(— + Ix^)) dx < a/tt/c exp(3/c). □ 
-87r/75 
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Proof of Theorem\^ Our approach will be to bound J \F{0, 0)| over a variety of regions 
whose union covers IZ^. To make the comparison of these bounds with J^/ F{6, 4>) easier, 
we note that 

I F{e,(f))ddd(f) = exp{0{m' + n'))lo = exp{0{m^' + n^'))h, (47) 
Jtz' 



where 



TT N 1/2 ^ / Nl/2 JL / vr \V2 



A„j n(4..) 
•• j=i k=i 

m/2 / 7J- \n/2 



^ / TT \^/^ / IT V 

^ ~ \AnJ \AmJ 



To see this, expand 



Ai, = An + a,, = An exp — ^-^ + ■ ■ 



2 

O" ■ 



and similarly for A.^, and compare the result to Theorem [2] using the assumptions of 
Theorem [1] It may help to recall the calculation following (H2|) . 

Take k = 7r/300 and define Xq, Xi, 0:299 by Xg = 2iK. For any i, let Si(i) be 
the set of {0, </>) such that 6j G [x^ — k,,X£ + k] for at least Km/ix values of j and 0^ ^ 
[— a;^ — 2/t, —xi + 2fi:] for at least n values of k. For (0, 0) G Oj + 0^ ^ ["'t? 't] for at 

least Kmn^/vr pairs (j, fc) so, by Lemma HJ (f))\ < exp(— CiAjnjnmn^) for some Ci > 

which is independent of C. 

Next define iS2(^) to be the set of {6, (j)) such that 9j G [x^ — k, x^ + k] for at least kui/ti 
values of j, 0^ G [—Xg — 2n, —Xg+2K\ for at least ra — ra^ values of and 6j ^ [x^ — Sk, x^+Sk] 
for at least imf values of j. By the same argument with the roles of 6 and reversed, 
\F{.^i 4>)\ ^ 6xp(— C2Ajnjnm^n) for some C2 > independent of £ when (0, 0) G iS2(£). 

Now define 7^i(£) to be the set of pairs {6, cf)) such that 6j G [x^ — 3/t, x^ + Sk] for at 
least m — values of j, and 0jt G [— — 3/t, — + Sk] for at least n — values of k. 
By the pigeonhole principle, for any 6 there is some I such that [x^ — k,,X( + k] contains 
at least Km/vr values of 6j. Therefore, 

299 299 

Since the total volume of ( IJ^^il^))'^ is at most (2m)™'^", we find that for some C3 > 0, 

/ \F{e,(j))\ded(t) 

< (27r)'"+"(exp(-C3Ajnin"^n'') +exp(-C3Aj^inm^n)) 

< e-"/i. (48) 



(u 
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We are left with {6, cf}) E IJf'^il^)- If subtract Xi from each 9j and add to each 
the integrand F{0, (p) is unchanged, so we can assume for convenience that i = and 
that {6,(l>) G TZi = 7?.i(0). The bounds we obtain on parts of the integral we seek to 
reject will be at least 1/300 of the total and thus be of the right order of magnitude. We 
will not mention this point again. 

For a given 6, partition {1,2,..., m} into sets Jg = Jq{0), Ji = Ji{6) and J2 = ^2(^)) 
containing the indices j such that \6j\ < 3k, 3k < \6j\ < 15k and \6j\ > 15k, respectively. 
Similarly partition {1,2, ...,n} into Kq = KQ{(f)), Ki = Ki{<p) and K2 = K2{(f)). The 
value of \F{6, can now be bounded using 

fjk{Oj + 

exp(-A^i„(% + + nAnin(^j + if (j, k) e (Jo U Ji) X {Ko U K,), 

< < Vl-4A^i„(l-cos(12«:)) < e-^-'"/^^ if (j, k) G ( Jq x K^) U ( J2 x K^), 

1 otherwise. 

Let 12(^2, be the contribution to J^JF^O, <p)\ of those {6,cf)) with IJ2I = "^2 and 
I -^2 1 = ^^2 • Recall that | Jg | > m — nf and | Kq \ > n — rf . We have 



/2K,^2)< (2vr) 

\m2j \n2j (49) 

X exp(--^A^in(n - n')m2 - ^A^^^{m - nf)n2)l2{m2, ^2), 

where 

/15k p15k 
■ ■ ■ I exp f-Ainin ^ + 0fc)^ + ^Ainin ^ + j dd' dcf)' , 

-15k ^ — 15k ^ J J / 

and the primes denote restriction to j G Jq U -^1 a-i^d k G Kg U . Write m! = m — m2 and 
n = n — n2 and define ^' = (m')~"^ 6*^, = 9j — O' for j G Jq U -^i? 0' = (^')~^ Y^k 
4>k = (j)k — 4*' for ^ £ -K'o U i^i , z/' = 0' + ^' and 5' = O' — 4) . Change variables from [6' , <p') 
to {9j I j G J3} U {0fc I k G -^"3} U {u, 6'}, where J3 is some subset of m' — 1 elements of 
Jq U Ji and X3 is some subset of n' — 1 elements of K^UKi. From Section H] we know that 
the Jacobian of this transformation is m'n /2. The integrand of l'2 can now be bounded 
using 

and 

Y!je, + < 27n' 6^ + 27m' + 27m'n'u'\ 
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The latter follows from the inequality {x + y + zY < 27(x^ + 1/^ + z'^) valid for all x, y, z. 
Therefore, 

0(1) r'^^'^ pSOk pSOk / ^ / V 



where g{z) = —z^ + l^;"^. Since g{z) < for |z| < 30k, and we only need an upper bound, 
we can restrict the summations in the integrand to j G J3 and k G K:^. The integral now 
separates into m' + n' — 1 one-dimensional integrals and Lemma E] (by monotonicity) gives 
that 



^(m +n )/2 



/^(m2,n2) = ^(^-W-i)/2, .„.V2-V ..V/2-i "^P(^^^V(-^--^^^+^V^-^-^°^^^^)- 
^min (m) / (n) ' 

Applying (jH]) and (jlO]), we find that 



m n 



E E ^2(^2,^2) =0(e-'=^^"' + e-'^^^")A (50) 



77x2=0 n2=0 



for some C4 > 0. 



We have now bounded contributions to the integral of |F(0, 0)1 from everywhere 
outside the region 

X ={{e,4)) \ l^jl, < 15/t for 1 < J < m, 1 < A; < n }. 

By Lemma m we have for (0, 0) G C'^^'^ (which includes A") that 

m n m n 

\F{d, 4>) I < exp (- 5^ ^ A,,(^, + 4 + z.)2 + i ^ ^ A,fc(e;. + 0, + i.)^) , 

j=l k=l j=l k=l 

where 6j = 6j — ^, 0^ = 0^ — and u = 6 + (f). As before, the integrand is independent 
oi 6 = — (j) and our notation will tend to ignore 6 for that reason; for our bounds it will 
suffice to remember that 6 has a bounded range. 

We proceed by exactly diagonalizing the (m+n+l)-dimensional quadratic form. Since 
YlT=i = Y2=i 'Pk = 0, we have 



m n m 



2 

.1^ 



j=l k=l j=l k=l 

m n m n 

+ 2 ^ ^ ajjj^k + 2z^ ^ ajjj + 2z/ ^ a.fc^fc- 

j=l k=l j=l k=l 
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This is almost diagonal, because = 0{n~^^'^), and we can correct it with the slight 
additional transformation (/ + X~^y)~^^^ described by Lemma [H where X is a diagonal 
matrix with diagonal entries Aj,, A,^ and A„. The matrix Y has zero diagonal and 
other entries of magnitude 0{n~^^'^) apart from the row and column indexed by u, which 
have entries of magnitude 0{n^^^). By the same argument as used in Section [4.21 all 
eigenvalues of X~^Y have magnitude 0(?t,~^''^), so the transformation is well-defined. 
The new variables i'&j}, {0k} ^-nd z> are related to the old by 



^1, • • • , ^m, 01, • • • , (Pn, l^f = {I + X ^Y) ^/^(t^i, . . . , (^1, . . . , (f^, 



We will keep the variable 5 as a variable of integration but, as noted before, our notation 
will generally ignore it. 

More explicitly, for some di, . . . , c?^, d'l, . . . ,d'n = 0{n^^^^), we have uniformly over 
j = 1, . . . ,171, k = 1, . . . ,n that 

m n 

q=\ k=l 

m n 

k = 0k + Y. + E O{n-^)0, + d{n-^'^)v, (51) 

j = l q=l 

m n 

V = i> + Y dj^j + ^ d'kcpk + 0(n~^)z>. 
j=i k=i 

Note that the expressions 0{ ) in flSTl) represent values that depend on m, n, s, t but not 
on {i^j},{(^fc},z>. 

The region of integration X is (m+n)-dimensional. In place of the variables {6, 4>) 
we can use {6, 0, z/, 6) by applying the identities 9^ = — Yl^=i 4>n = ~ Ylk=i k- 

(Recall that and (j) don't include and The additional transformation (15T!) 

maps the two just-mentioned identities into identities that define and 0^ terms of 
(■j9, (p, z>), where "d = {^i, . . . , ^m-i) and = {0i, . . . , These have the form 

m— 1 71—1 

, (52) 

m— 1 n— 1 ^ ' 

j=i k=i 

Therefore, we can now integrate over (i9, z>, 5). The Jacobian of the transformation 
from {6,4)) to {0,(f),i',6) is mn/2, as in Section HI The Jacobian of the transformation 
T4^{'&, (p, v) = {6, 4>, v) defined by fl^ can be seen to be 1 + 0{rr^^'^) by Lemma [21 using 
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the fact that the oo-norm of the matrix of partial derivatives is 0{n ^^'^). This matrix 
has order m + n — 1 and can be obtained by substituting (152]) into (15T]) . 

The transformation T4 changes the region of integration only by a factor 1 + 0{n~^^'^) 
in each direction, since the inverse of (15T|) has exactly the same form except that the 
constants {dj}, {d'j^}, while still of magnitude 0(n~^^^), may be different. Therefore, the 
image of region X lies inside the region 

y = {{'d, (p, v) I \^k\ < 31k (1 < j < m, 1 < A; < n), |z>| < 31^: }. 



We next bound the value of the integrand in 3^. By repeated application of the 

1 2 _|_ 1 



inequality xy < + ^y^, we find that 



m n m n 

j=l k=l j=l k=l 



where we have chosen | as a convenient value greater than |. Now define h{z) = —z^ + ^z'^. 
Then, for (p, z>) G 3^, 

m n 

\F{e,(f>)\ < exp(^A,.M^,) + + A../i(z>)) 

j=i k=i 

m—l n—1 

< exp( ^ A^,h0j) + Y,A.kh{^k) + A.,h{u)^ (53) 
j=i k=i 

m—l n—1 

= exp{A,.h{u)) Yl exp{Aj,h0j)) JJ exp(A.fc/i(^,,)) , (54) 
j=i k=i 

where the second line holds because h{z) < for \z\ < 31k. 
Define 

Wo = {{'d,p,u)ey\ < K'^'^' (1 < J < m - 1), 

\0k\ < (1 < A; < n - 1), 

m—l n—1 

j=l k=l 

Also define similar regions Wq, W'l, W'2 by omitting the variables di, ipi instead of d^, ipn 
starting at (153|) . Using (ISTI) . we see that T4, and the corresponding transformation that 
omits i^i and (fi, map 7^ to a superset of Wq H W2 fl W( fl W2. Therefore, A* — 7^ is 
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mapped to a subset of Wi U (Wq - W2) UW[U (Wq - W2) and it will suffice to find a 
tight bound on the integral in each of the four latter regions. 

Denoting the right side of by Fq^-O, (p, z>), Lemma [5] gives 

/ FQ{'d,<^,u) d'dd(pdu = exp{0{m' + n''))lQ. (55) 
Jy 



Also note that 

r-3lK 



exp{ch{z)) = 0{l)exp{ch{zo)) (56) 
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for c,Zq>0 and Zq = o(l), since /i(2;) < H^Zq) for -Zq < ;z < 31k. By applying (156|) to each 
of the factors of (IMl) in turn, 

[ Foi^, (p, z>) dm(pdu = O (e-^^^'"'' + e-'^^^"'') Jo (57) 

for some Cg > and so, by (1^ and (IFTI) . 

/ Fo(i9, (p, z>) dM(pdu = exp(0(m^ + n^))/o. 

Applying Lemma [3] twice, once to the variables di-d^, . . . , (i^_i'i?m„i, (i'i<^i, . . . , d'n-i0n-i 
and once to their negatives, using M = 0{n~'^), N = m + n — 2 and t = n~^^^, we find 
that 

/ Fo{'d,<p,u)d'dd(pdu = 0{e''''^^) [ Fq{^,(p,v) dMCpdv 

= 0(e-"''Vo. (58) 

Finally, parallel computations give the same bounds on the integrals over VV( and 
Wo - W2. 

We have now bounded j\F{0, (p) \ in regions that together cover the complement of TZ. 
Collecting these bounds from fHHj) . fl50|) . fl57|) . fl58|) . and the above-mentioned analogues 
of fl571) and fl58|) . we conclude that 



/ \F{e, (f))\d9d(j) = 0{e-''^'^'^ + e-'^^"'')/o 
for some Cy > 0, which implies the theorem by ( H7I) . □ 
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Appendix: Estimating an integral 



In this appendix we estimate the value of a certain mult i- dimensional integral. A similar 
integral appeared in [7] and variations of it appeared in [H El [6]. However, none of 
the previously published variations meet our present requirements entirely. We will meet 
them here, and also introduce a new method of proof that gives a better error term. 

It is intended that this appendix be notationally independent of the rest of the paper. 
We have used new symbols where possible, but even in the few remaining exceptions, 
assumptions about the values of variables stated earlier do not apply here. 

Theorem 4. Let e' , e" ,e"' ,e, A be constants such that < e' < e" < e" , £ > 0, and 
< Z\ < 1 . The following is true if e" and e are sufficiently small. 

Let A = A{N) be a real-valued function such that A{N) = Q{N''^ ). Let dj = dj{N), 
Bj = Bj{N), Cjk = Cjk{N), Ej = Ej{N), Fj^ = Fj^iN) and Jj = Jj{N) be complex- 
valued functions (1 < j,k < N) such that Bj, Cj/., Ej, Fjj^ = 0{N^), dj = 0{N^^'^'^^), and 
Jj = 0{N~^^'^^^), uniformly over 1 < j, k < N . Suppose that 



. N N N N 

f{z) = exp f -iAT Y.^' + Y. ^^oA + ^ E + E ^^^^3 

N N N V 

+ ^ E ^^-^J + E ^^-^^^^ + E '^I'l + ^(^) ) 

j=l j,k=l j=l 

is integrable for z = {zi, Z2, ■ ■ ■ , Zj^) G and 5{N) = max^^zu^lS^z)] = o(l), where 
Un = {z CR^ \ \zj\ < AT-i/^+e l<j<N}, 

where i = i{N) satisfies e" < 2e < e" . Then, provided the 0() term in the following 
converges to zero, 

/ f{z)dz={^] exp{0,+ 0, + O{{N-^ + 6{N))Z)), 
Jun \AJ\I / 



2 

Zk 



where 



1 ^ 1 ^ 15 ^ 3 ^ 

&i = — - — > dj H — ^ y dj H ^ — > Bj H ^ — - y BjCjk 

nAAT^-^ ■' Aa2t\t2/_^ J a A-i AT ■' Q A'i AT^ ^—^ ■' ■' 



2AN ^ ' AA^N^ ^ ' WA^N ^ ' SA-'N . ^ , 
j=i j=i j=i j,k=i 



1 ^ 2 Af 1 ^ 



j,k,i=l j=l J,k=l 
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O — ^ + a E + \^ a 



1 ^ . 3 ^ . , 1 ^ , , 

+ W¥ ^^^^ + "'^^^^■^ + ii^ ^ ^^-^^ + 

1 ^ . . 3 ^ 



exp 



^ .7=1 :;=i 

N 1 ^ \ 



Proof. Our method of proof will be integration over one variable at a time. This method 
is conceptually simple but technically challenging. Assistance from a computer-algebra 
system is recommended. 

Let Hj^ j^ j^ be a functions of N for each 1 < ji,j2, ■ ■ ■ ,jk ^ ^ l^t Pi,P2, ■ ■ ■ ,Pk 
be non-negative integers. Let 1 < j < A^ + 1. Define the generalized moment 

jl,--;jk 

where the summation is over 

{(Ji, ■■■,jk)\ ■ ■ -JkM = A;, 1 < J, < J - 1 if = 0, J <ji<N Hp, > 0}. 

We will customarily omit the argument z as it will be clear from the context. Note that 
the indices ji, ■ ■ ■ , jk are reserved to this notation and always index the position their 
name suggests; for example 

^3 



r/,(a,j3,0) = E 



i<i2<i-i 

We also need the defective moment Vj{Bj_^^ jJ P11P21 ■ ■ ■ iPk) which is the same as 
rij{Hj^ jj Pi,P2, ■ ■ ■ ,Pk) except that the index value j is forbidden; that is, the condition 
ji, ■ ■ ■ ,jk 7^ j is added to the domain of summation. 

Some properties of these moments that we require are listed below. Assume that 
z eUj^. Then 

\vJiH,,,...Jp^,...,Pk)\<ras^\H,^,...J n'^^-'/'^'^^^^-^^'^\ (59) 

?7i(-f^ji,...jfc I Pi, ■ ■ ■ ,Pfe) = if Pi = for any i, (60) 

?7jv+i(^ii,...jfe I Pi, ■ ■ ■ ,Pfe) = if pj > for any i, (61) 
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i\Pi>0 

(62) 

Vj+i(Hj^,...,jJPi, ■■■,Pk)^ v'j{Hj^,...,jJPi, ■■■,Pk) 

i\Pi=0 

The last two equalities require j < N. 

The product of generalized moments ^7j(-Pjj,...j,,| Pi, • • • ,Pk) and rij{Qj_^^ j^ \ gi, . . . , Qi) 
can be written as a sum of generalized moments. Define $ to be the set of injec- 
tions : {1, 2, . . . , £} ^ {1, 2, . . . , + £} such that (a) 0({1, 2, . . . , £}) U {1, 2, . . . , fc} = 
{1, 2, . . . , 101} for some integer |0| depending on 0, (b) ior 1 < i < j < i, if 0(i), > k 
then (f){i) < (t){j), and (c) for 1 < i < £, q'j = (0(i) > A; or = 0). For e $ and 
1 < « < |0|, define rj = Pj + where the first term is omitted iii > k and the second 

term is omitted if i is not in the range of 0. Then 

Pi' ■■■^Pk) Vj{Qh,...,je\ • • • ' 

For example, 

^i(«iiJ2l 0' 2) r;,(/3,-^ I 3) = Vjiaj^jAi I 0' 2, 3) + ^.(".ij/i^ I 0' 5), 

where the two terms correspond to the injections 0(1) = 3 and 0(1) — 2. Exactly the 
same formula holds for defective moments. 

For 1 < j < AT + 1, define 

F^{z) = v,{-AN + a,J 2) + v^{B^^N + ^^-J 3) + r/,(4iV + 4,.J 4) + v^iJ^J 1) 
+ ^,(4,^ I 2, 1) + ^,(4,^ I 2, 2) + r/,(ro| 0) + ^,(ro,o| 0, 0) + 77,(^,0! 1, 0) 

+ ^,(r2,ol 2, 0) + ^,(ro,o,ol 0, 0, 0) + ^,(A,o,ol i, o, o) + Vjir,^i,o\ i, i, o) 

+ r^,(r2,o,o| 2, 0, 0) + Vj{r2,i,o\ 2, 1, 0) + Vj{r2,2,o\ 2, 2, 0) + r^,(A,o,o,ol 1, 0, 0, 0) 
+ ^,(A,i,i,ol 1, 1, 1, 0) + ?7,(^2,i,o,ol 2, 1, 0, 0) + ?7,(^2,2,i,ol 2, 2, 1, 0), 
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where 



^1,1,1,' 







J2 



° ~ 2AN AA^N"^ QA^N^ AA^N AA^N IQA^N 

°'° 4i^Ar2 4i2^2 g^3^2 ^^3^3 g^4^3 

2iiV 2i2iV^ lei^iV^ 2i^iV''^ 2i-'^A^^ ' 

J- 2 — ^^^T 1 ~ 1 1 ^^^o T 1 o 1 

AA^N 2AN 2AN 2A^N^ 2A^N^ 

p _ ^iii-ihh _|_ "'is ilia i2J3 

P _ i^h + ^i3)^i3i2^i3ii I '^iaii '-^i3i2 

7-1 _ "^^i2i3 ^iaii -^ia I (-^iiia ^i3h)^32h 
2,1,0 ^ 2i3;v2 + 2i2iV^ 

n n h C C 
P _ ^iaii i3i2 I "'ia i3ii iai2 

P _ ^hji^uh^uh I ^i4ia ^J4i2 ^iii2 

1,0,0,0 - 14 Ar4 ^ o 14, r4 ' 



AAA 



A A A A A 

P _ ^i2i4^i4ii^i4i3 , ^hh^iih hh 
2,1,0,0 ~ TTsTTs ' TTws ' 

AAA 

7-1 _ ^i3i4 i4ii i4i2 
4i^7V2 ■ 

Note that FAz) is independent of for i < j. The key properties of Fj{z) for z E 
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are 



/(z) = exp(Fi(z) + 5(z)), (65) 
/ . exp(F,.(2)) dz^ = exp{F^^,{z) + 0(iV-i-^)) (j < N). (66) 

Equation 0651) is easily seen after applying (H7]) to eliminate most of the terms. Proof 
of dSSD requires a tedious calculation which we now outline. 

First, apply (l62l) to make explicit the dependence of Fj{z) on Zj (a polynomial of 
degree 4). Then expand 

exp(F,.(2)) = exp(i?o(^)) exp(-iArz|) [l + R^{z)z^+ R^{z)z] + ■ ■ ■ + 0{N-^-^)) , (67) 

where each Ri{z) is independent of Zj and contains defective moments only. As seen by 
applying fl59p . only a finite number of terms are required to achieve the requested error 
term. The factor exp{—ANz'^) comes from the first term rjjl—ANl 2) of Fj{z). Products 
of moments that occur need to be rewritten as sums using (IMl) . 

Next, integrate (1671) over Zj using 

/ /V^^^'dz= — S^-^,/^(l + 0(exp(-ciV^^"-'))), 

for fixed A; > 0, for some c > 0. Here we have used the assumptions that A = f2(A^~^ ) 
and e' < e" < 2e. The result of the integration has the form 



exp(i?o(^)) (1 + S{z)) = exp{R,{z) + log(l + S{z))). 



AN ^ ^ - ^ ^ - Y^jV 

Since S{z) = o(l) (in fact S{z) = 0{N^'^^'^^^^'^ ^'^^) for some A;), we can expand the loga- 
rithm using fl59|) again to limit the expansion to finitely many terms. Finally, apply fl63l) 
to rewrite the defective moments in terms of ordinary generalized moments. The result 
is the right side of ( !66|) . 

If all the coefficients dj, Bj,Cjf., Ej, Fji^, Jj were real, we could apply fl65l) and fl66l) 
immediately to find that 

f{z) dz = l^-^j exp(F^+i + 0{6{N) + iV'^)), (68) 

noting that Fj^_^_i{z) is independent of z. 
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When the coefficients are complex, we must take more care. Equation (165|) only allows 
us to write 



f{z)dz= [ exp{F^{z))dz + 0{5{N)) [ \exp{F^{z))\ dz. 



(69) 



Let F*{z) be the same as Fj{z) except that the coefficients a^, Bj, Cj^, Ej, Fj)., Jj are 
all replaced by their real parts. Clearly 

|exp(Fi(2;))| =exp(Fi*(2)), 

and so, as in fIBSl) . 

\N/2 



exp(Fi(z)) d2 = 0(1) 



TT 



AN, 



exp(Fj 



N+l) 



(70) 



From fl66l) we have for 1 < j < that 

-1/2 

|exp(F^+i(2)) 



TT 



< 



AN 



TT 



AN 



(i + o(iv-^-^)) 



exp(Fj(z)) (izj 



-1/2 



(1 + OiN-^-^)) j |exp(F,(z))| 



so we have by induction starting with (|70|1 that 
y |exp(F_,+i(2))|rfZj+i---dz^ < 

Returning to (!66|l . we find that, for 1 < j < A^, 
y exp(Fj(2)) (i^;^ ■ ■ ■ dzj~^ = 



exp(F^+i). 



AN 



7r 



AN 



+ 0{N 



j exp(F^.+i(;2) + 0(Ar"^-^)) rfz^.+i ■ ■ ■ 

y exp(Fj+i(z)) ■ • • d^^v 

j I exp(Fj_|_i(2)) I ■ ■ ■ dz 



TT 



y exp{Fj_^i{z)) dzj+i ■■■dzN 

^ x(JV-j+l)/2 



AN 

+ OiN-^-^) 



AN 



) 



exp(F;^+i). 



By induction on j, this gives 



/ exp(Fi(z)) dz=[^\ (exp(F^+0 + 0(iV-^) exp(n+0), 

\AN j 
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which, together with (!69|) and (ITOl) gives 

/ f{z) dz=(^] (exp(F^+0 + 0(iV-^ + 5{N)) exp(n+i)) 

= (^^j exp(F^+i + 0(iV-^ + 5(iV))Z), (71) 

where Z = exp(F^_,_i — Re(F/v+]^)) and the last hne is vahd if (A^~^ + 6{N))Z = o(l). 
Applying (l6Ti) to the definition of Fj{z), we find that 

FN+i=O, + 0, + OiN-^), 

from which it follows that Z has the value in the theorem statement to within a multiplied 
constant. Also note that Z > 1, which is easiest to see by noting that the argument of 
the exponential is a non- negative quadratic form for each j. The theorem now follows 
from ([7ID. □ 
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